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ABSTRACT 

We present a direct detection of the growth of large-scale structure, using weak grav- 
itational lensing and photometric redshift data from the COMBO-17 survey. We use 
deep i?-band imaging of two 0.5 x 0.5 square degree fields, affording shear estimates 
for over 52000 galaxies; we combine these with photometric redshift estimates from 
our 17 band survey, in order to obtain a 3-D shear field. We find theoretical models 
for evolving matter power spectra and correlation functions, and fit the correspond- 
ing shear correlation functions to the data as a function of redshift. We detect the 
evolution of the power at the 7.7<r level given minimal priors, and measure the rate 
of evolution for < z < 1. We also fit correlation functions to our 3-D data as a 
function of cosmological parameters <7s and Ha- We find joint constraints on fl\ and 
us, demonstrating an improvement in accuracy by a factor of 2 over that available 
from 2D weak lensing for the same area. 

Key words: gravitational lensing — cosmology: cosmological parameters, dark mat- 
ter, large-scale structure of Universe 



1 INTRODUCTION 

The evolution of the distribution of matter is a topic of cen- 
tral concern to cosmology. We seek to understand the history 
of matter environments at all epochs, from the near uniform 
distribution of matter at z ~ 1000, to the highly skewed dis- 
tribution at present. Detailed knowledge of this evolution 
will lead to precise measurements of cosmological parame- 
ters such as the dark energy density fi„, and its equation of 
state w (e.g. Refregier et al 2003, Benabed & van Waerbeke 
2003). Furthermore, a comparison of the evolution of dark 
matter and baryonic matter will afford a deep understanding 
of the formation of baryonic structures. 

There are several approaches available for studying 
matter evolution. For example, one can compare the matter 
power spectrum inferred from the Cosmic Microwave Back- 
ground with that inferred from low-redshift galaxy surveys 
(e.g. Percival et al 2002, Spergel et al 2003); or if one is 
principally interested in baryonic structure evolution, one 
can survey galaxy counts out to higher redshifts and directly 
observe evolution in the count power spectra (eg Phleps & 
Meisenheimer 2003). Here, we will study another approach: 



direct measurements of matter power spectrum evolution 
via weak gravitational lensing combined with photometric 
redshifts. 

Weak gravitational lensing has proved to be a highly 
useful method of obtaining cosmological information. The 
weak lensing phenomenon is observed as the slight align- 
ment of neighbouring galaxy images, due to light rays being 
deflected by gravitational potential fluctuations along their 
light paths. As the deflection is due to the gravitational en- 
vironment through which the light passes, the phenomenon 
is sensitive to all matter, both dark and baryonic. Lensing 
therefore allows measurement of the distribution of the en- 
tire matter content of the universe. 

The use of weak lensing (often without redshift esti- 
mates for individual galaxies) to constrain the matter power 
spectrum is now well developed. The shear estimates for 
large numbers of galaxies have been used to measure the 
statistics of an effective 2-dimensional shear field; theoret- 
ical calculations of these statistics can be fit to the data 
in order to estimate cosmological parameters (see e.g. van 
Waerbeke et al 2001, Hoekstra et al 2002, Bacon et al 2003, 
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Brown et al 2003, Jarvis et al 2003) or the statistics can be 
used to measure the dark matter power spectrum (e.g. Pen 
et al 2003, who measure the 3-D power spectrum at a fixed 
redshift). 

Now, with the increase in availability of reliable photo- 
metric redshifts, there is a great deal of interest in combin- 
ing redshifts with weak lensing. For example, Wittman et al 
(2001, 2003) have used shear and redshifts to measure the 
mass and redshift of two clusters directly from the 3-D shear 
field. In addition, several recent theoretical studies examine 
how to directly reconstruct the full 3-D gravitational poten- 
tial from lensing together with photometric redshifts (Taylor 
2001; Hu & Keeton 2003; Bacon & Taylor 2003); this will 
allow examination of peak statistics to obtain cosmological 
parameters, and studies of galaxy formation as a function of 
environment. 

A more directly statistical approach to the shear esti- 
mators is to use redshifts to permit shear power spectrum 
tomography (e.g. Seljak 1998, Hu 1999, 2002, Huterer 2002, 
King & Schneider 2002b, Heavens 2003). In this methodol- 
ogy, the galaxies are divided into several redshift bins, in 
each of which shear correlation functions or power spec- 
tra are measured; by comparing these to theoretical mod- 
els, cosmological parameters can be estimated. This paper 
will utilise an extension of this approach. This methodology 
also forms the starting point for removal of systematic pollu- 
tants of the cosmological shear signal such as intrinsic galaxy 
alignments (Heymans et al 2004, Heymans & Heavens 2003, 
King & Schneider 2002a,b). An alternative 3-D statistical 
method (particularly for measuring w and to') is provided 
by Jain & Taylor (2003), which is further investigated by 
Bernstein & Jain (2004). 

In this paper, we seek to measure directly the evolu- 
tion of the matter power spectrum. We will achieve this 
using a maximum likelihood fit between data and theoret- 
ical power spectra evolving with redshift. We will also fit 
evolving power spectra expected for various values of cos- 
mological parameters, in order to demonstrate the power of 
this method for precise determinations of dark energy pa- 
rameters with future large lensing surveys. 

We will use the COMBO-17 survey (Wolf et al 2001) to 
measure the dark matter evolution. This survey has already 
been extensively studied from a weak shear perspective; in 
particular, it has yielded precise measurements of mass-to- 
light correlations for a supercluster (Gray et al 2002), mea- 
surements of the shear power spectrum arising from large- 
scale structure (Brown et al 2003), and 3-D gravitational 
potential maps (Taylor et al 2004). One of the great assets 
of the COMBO-17 survey is the existence of accurate pho- 
tometric redshifts (Az ~ 0.05) for galaxies with z ^ 1. Our 
previous statistical lensing study with COMBO-17 only used 
these redshifts to obtain median redshifts for lens and source 
planes; in this study, we will use the redshifts for each indi- 
vidual galaxy to compare each galaxy pair's distortion with 
what would be expected theoretically for galaxies at these 
3-D positions. 

The paper is organised as follows. In section 2, we will 
describe the construction of the theoretical power spectra 
which we require. In particular, we discuss the modeling of 
the underlying evolving matter power spectrum; we go on to 
discuss the shear power spectra expected given this matter 
power spectrum. This includes not only power spectra for 



shear in a redshift shell, but also for cross power spectra for 
shears at two different redshifts. We also find inversion equa- 
tions for moving between shear and matter power spectra, 
and construct phenomenological models to measure evolu- 
tion of the power spectra. 

We proceed in section 3 to describe the COMBO-17 sur- 
vey, giving details regarding the observations, photometric 
redshift estimation and weak shear analysis. 

Section 4 describes the analysis of the data. We measure 
the rate of evolution of the matter power spectrum using a 
simple model. We also find constraints upon the cosmologi- 
cal constant and amplitude of the power spectrum by fitting 
models with varying values of these parameters to the 3-D 
shear field data. 

Finally, in section 5 we discuss the implications of our 
results, and discuss the future development of the method 
on large shear datasets. 



2 POWER SPECTRA AND CORRELATION 
FUNCTIONS 

In this section we will describe our approach to calculat- 
ing theoretical models for the matter power spectrum and 
shear power spectrum for various cosmologies. We will use 
these models to compare with the data, in order to mea- 
sure cosmological parameters. We also describe here how we 
obtain cross-correlation models for shear estimators at dif- 
ferent redshifts, and how one can calculate the matter power 
spectrum from the shear power spectrum and vice versa. Fi- 
nally, we will present a phenomenological model suited to a 
direct measurement of the evolution of the power spectrum. 



2.1 Matter Power Spectrum 

We begin by describing the method used for calculating the- 
oretical matter power spectra for different values of cos- 
mological parameters. This will require us to calculate a 
full non-linear power spectrum, which we will obtain via 
the intermediate step of calculating the linear power spec- 
trum. First then, we require a model for the linearly evolv- 
ing power spectrum A 2 (A:) in logarithmic increments of co- 
moving wavenumber k. In order to achieve this, we use an 
initially Harrison- Zeldovich power spectrum, evolved with a 
transfer function T, growth factor g described below, and 
a (1 + z)~ 2 factor as expected for Einstein-de Sitter linear 
growth of a matter-dominated expansion (c.f. Bardeen et al 
1986): 



A 2 (M) 



Ak 3 T 2 (k)k 



(l + z)V(0)' 



(1) 



where A is a numerical constant which we will normalise 
presently. The transfer function T for the Cold Dark Matter 
scenario is given by (e.g. Bardeen et al 1986; c.f. Efstathiou, 
Bond & White 1992, Bond & Efstathiou 1984): 



In (1 + 2.34?) 
TGDM(fc) = 2.34 g X 



(2) 



[l + 3.89q + (16.1q) 2 + (5.46g) 3 + (6.71g) 4 ] 



41 -1/4 
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with q = k/hT. Here, h is the present-day Hubble constant 
Ho divided by 100 km s _1 Mpc -1 , and the 'shape parameter' 
r is defined as (Sugiyama 1995) 

r = fi m ,oftexp(-n B ,o(l-v / 2ft/fi m ,o)). (3) 

This equation introduces the present-day density parameters 
for matter, Q m ,o, and baryons, £Ib,o- In order to calculate 
A 2 , we define a growth parameter g which describes the 
difference in growth between a particular cosmology and an 
Einstein-de Sitter cosmology at a given redshift z (Carroll 
et al 1992): 



(i + JLn„(z)) 



(4) 
(5) 



In turn this equation requires the evolution equations for 
the density parameters of matter, fi m , and of the vacuum 
energy, 



fim(«) 



o + (l 



a)£l m ,o + (a 3 



(6) 



(7) 



a + (1 — a)Q m ,o + (a 3 — a)f2„,o ' 

where the expansion parameter a = 1/(1 + 2). Now we have 
the equipment for calculating the linear power spectrum, we 
can progress to estimating the non-linear evolution of the 
power spectrum, via the formalism of Smith et al (2003). 

This method calculates the fully non-linear power us- 
ing the insights afforded by the Halo Model (e.g. Ma & Fry 
2000, Peacock & Smith 2000, Seljak 2000). The power is de- 
composed into two terms, a quasi-linear term describing the 
power from the clustering of halos, and a halo term describ- 
ing power from the clustering of objects within each halo. 
Each of these terms is calculated by Gaussian filtering the 
linear power to find a non- linear threshold wavenumber k a , 
spectral index and spectral curvature. These are used to find 
11 coefficients which are functions of the calculated spectral 
properties; the quasi-linear and halo power terms are then 
found as functions of the linear power and these coefficients, 
finally resulting in the full power A 2 . Full details of this 
procedure can be found in Smith et al (2003). 

Finally we move from the dimensionless power A 2 to 
the standard power spectrum P via 



P NL (k, z) = 2n 2 A 2 NL (k, z)/k 3 . 



(8) 



We normalise the power spectrum at the present epoch 
by setting a value for the quantity erg, measuring the rms 
fluctuation of density when smoothed by a top-hat filter 
with radius R = 8/i _1 Mpc: 



2 

a s = 



k 2 dkP(k,0)- 



:(sin(kR) - kRcos(kR)). (9) 



(kR) 3 1 

Thus, for particular choices of as, f^m.o, ^u,o, h, Qf, and 
fife , we have a model for the evolving non-linear power spec- 
trum at any required redshift. 



2.2 Shear Power Spectrum 

Now that we have a model for the matter power spectrum at 
any epoch, we can proceed to calculate the 3-D shear power 



spectrum and correlation function. It is convenient to obtain 
these via the effective convergence k, a measure of the sur- 
face mass density, which is given by (following Bartelmann 
& Schneider 2001, Ch 6 throughout this section): 



.„ 3Hjn m p S[f K (r)e,r] 
= o„7 / drW{r)f K (r) 



2c 2 



a(r) 



(10) 



Here, 9 represents the position on an image, 8 is the density 
perturbation field, c is the speed of light, m is the comoving 
radial distance to the horizon, and r is the comoving radial 
distance to a given mass concentration: 

r ( 2 ) = TT / [an m (0)+a 2 (l-fl m (0)-fl v (0))+a 4 Q v (0)r 1/2 da.(ll) 

In the equation for k above we also require the radial 
coordinate distance fx, 

( K~ 1/2 sin(K 1/2 r) (K > 0) 

f K (r) = { r (K = 0) (12) 

[ (-A:)- 1/2 sinh[(-^) 1/2 r] (K<0) 

with the curvature K given by 

K= (J^j (n m + 0,-l). (13) 

We also require the weighting factor W , 

W(r)= fK f (ra ~ r) H(r a -r), (14) 

where H is the Heaviside function and r s is the comoving 
distance to a particular source galaxy. 

We can now find the cross-power between convergence 
at any two redshifts. In order to do this we use Lim- 
ber's equation, which can be formulated (e.g. Bartelmann 
& Schneider 2001) to state that if we have two projections 
Pi,i = 1,2 of the density field S, such that 



/ 



Pi (0) = / dr qi (r)6[f K (r)8,r] 



(15) 



for some function qi(r), then the cross power spectrum will 
be 



fMr) 



fK(r)- 



(16) 



If we wish to calculate the cross power of k between two 
redshifts, then we see by comparing equations (10) and (15) 
that we do have two projections of the necessary form, with 
9i given by 



M _3# 2 o fK(r)f K (n-r)H(n-r) 
U ' ~ 2 c 2 " m f K ( n ) a(r) : 



(17) 



so we therefore obtain (noting that the two-point statistics 
of k and 7 agree, e.g. Blandford et al 1991), 



Pjl2(£, zi, 22) 



9#o » 2 

4c 4 



Jo 



dr 



- r) f K (r 2 - r) 
/jf(ri) //f(r 2 ) 



-Ps 



a 2 {r) \/if(r) 



(18) 



where I is the angular wavenumber. In the case where we 
wish to examine the power of the shear at one particular 
redshift, we can simplify this to 
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4c 4 



Jo 



dr 



f%{n 



Ps 



f 2 AnW(r) °\f K (r) 



.r .(19) 



We can obtain corresponding cross-correlation functions 
for shears at two different redshifts zi , z 2 . We will use three 
different correlation functions: Ci represents (7171), where 
7™ ' represents the first shear component of two galaxies, in 
a coordinate frame where zero position angle lies along the 
line joining the galaxies (e.g. Bacon et al 2003). C 2 repre- 
sents (7272)1 an d C — Ci + C2. Given these definitions, the 
correlation functions are simple transforms of the power: 



£d£ 
2ir 



P- l i2(t,zi,z 2 )Jo(£0), 



(20) 



Ci(0,zl,z2) 



47T 



47T 



P 7 i2 (£, z u z 2 )[J {16) + J 4 (16)] , (21) 



P-yi 2 (£,z u z 2 )[Jo(£9) - Jt(tO)]. (22) 



C(6,zl,z2)= 

Jo 

-f 

Jo 

C 2 (6,zl,z2) = ( 
Jo 

These are the 3-D shear correlation functions we have 
been seeking; we will use these functions to compare our 
data with various evolving cosmological models. 

2.3 Power Spectrum Inversion 

It is convenient to have a simple means of calculating the 
shear power spectrum from the matter power spectrum; this 
has been described in the last section. However, it is also use- 
ful to have a means of calculating the matter power spec- 
trum given the shear power spectrum; here we outline how 
this can be achieved for theoretical models. Noisy data are 
difficult to invert with this procedure; it would be more con- 
venient to fit models to the data and then use the results of 
this section to invert these models. 

In order to show how to invert the shear power spec- 
trum, we start by considering an integral of the form 



A(r)= / dr'B(r',r), 
Jo 



(23) 



where B(r',r) is a smooth, continuous function. Carrying 
out a partial differentiation of A with respect to r, we find 



dAjr) 
dr 



f 

Jo 



or 



(24) 



This result will be used below in order to untangle the mat- 
ter power spectrum from its integral projection found in 
calculating the shear. We also require the result 



d 
dri 



,Mn - r) 



Mr) 



(25) 



which can be easily verified directly from equation (12). Us- 
ing these two results repeatedly upon the shear power spec- 
trum P-y(l, r), we find that we can calculate the matter power 
spectrum: 



Ps(k,n) = 



4c 4 



a 2 (n) 



if 



9H*Q 2 m 2f K (n) \dr 2 



+ K 



f 3 ( r ■,dPj(f K {ri)k,ri) 
JK(ri) — 



dri 



(26) 



Alternatively, we can find a similar means of calculating 



the matter power spectrum from the cross power spectrum 
of shear at two different redshifts: 



Ps(k,n) = 



4c 4 



a 2 (ri)fK{r 2 ) 



d 



QHfol fK(n)f K (r 2 - n) dri 
dP 1 i 2 (f K (ri)k,r 1 ,r 2 ) 



fHn)- 



dri 



(27) 



These equations therefore permit us to find the matter 
power spectrum given a model for the shear power spec- 
trum. Note that upon noisy data, the first and second dif- 
ferentials in this equation can lead to unphysical negative 
power spectra. Thus we reiterate that a better approach is 
to fit a continuous shear model to the data which can then 
be inverted with these equations. 

2.4 Power Spectrum Model for Slope 
Phenomenology 

We conclude this section with a discussion of models for the 
matter and shear power spectra which allow us to directly 
examine the redshift evolution of the power spectrum. 

Firstly, suppose we have measured the shear cross- 
correlation function C between many redshifts, from which 
we can calculate the shear power spectrum. We will initially 
restrict ourselves to power within a redshift shell, P 7 ; we will 
consider P 7 i2 later. We will attempt to use our formalism 
with a simple model for the shear power spectrum: a power 
law in both the angular and redshift directions, i.e. 



P^{l,r) = Atr 13 = Ak a r a+f) . 



(28) 



The correlation function corresponding to this model 
calculated from equation (20) is: 

A2 a - 1 aT{a/2) a _ 2 _ a g 



(7(0, r) 



7iT(-a/2) 



(29) 



(One should note more generally that, if P 7 oc k a f(r), then 
C ot9- 2 - a r- a f(r).) 

Hence, if we fit a power law to the correlation function 
in the angular and redshift directions, we can immediately 
obtain estimates of the corresponding power-law approxi- 
mation to the shear power spectrum, and an approximation 
to the underlying matter power spectrum. To achieve the 
latter, we use the inversion equation (26); we find that the 
matter power spectrum is given (in a flat universe) by 



4c 



«„2/-_i r ,a+/3-l 



Ps(k,r) = - m ^A{a+(3){a+(3+\){a+(3+2)k a a z {r) 

This is a very useful formula; if we can adequately fit 
the shear power spectrum as a power law in the angular and 
redshift directions, we can instantly infer the corresponding 
power law matter power spectrum. 

We can attempt a similar approach with the cross power 
spectrum, P 7 i2, with a power law model keeping one redshift 
fixed: 



P-ri 2 (l,r u r 2 ) = Arr? = Ak a r° 



(31) 



However, this is less successful, as we obtain the follow 
ing form for the matter power spectrum: 

4 2/ \ a+P—1 

t a (ri)r 1 



Ps{k,n) 



4c 4 



■A(a+f3)(a+f3+l)r 2 k° 



r 2 - n 



(32) 
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This model fails to be physical; besides an unsightly 
divergence at = ri, we find that the matter power at 
the present epoch is either zero or infinite, depending on 
the power law slope. We therefore conclude that a power 
law model for the cross shear power is unsuitable for di- 
rectly calculating the matter power spectrum. Equivalently, 
we should note that the underlying matter power spectrum 
which produces a power-law cross shear power spectrum is 
unphysical. 

In order to find a measure of matter power spectrum 
evolution, able to take into account all the information from 
the shear cross power spectrum, one can take an alterna- 
tive approach. We will assume a very simple form for the 
underlying matter power spectrum, 

Ps{k ' ri) ^M§^ Akae ~ SZ ' (33) 

and will restrict ourselves to an underlying ACDM geometry 
in moving from this matter power to the shear predictions. 
This allows us to directly examine the matter evolution, 
rather than examining both matter and geometry change. 
The form chosen for Ps has the advantage of being positive 
definite as required for power; it is also a good fit to ACDM 
evolution in the non-linear regime that we wish to study. 

In order to remove some of the freedom from the above 
equation, we note that a — —1.2 provides an excellent fit 
to the angular dependence of the ACDM shear correlation 
function throughout z < 1 and 1' < 6 < 30'. We therefore 
allow only the matter amplitude A and redshift evolution 
slope s to vary. We will fit the corresponding shear corre- 
lation function calculated from equation (18) to the data, 
allowing us to directly measure the rate of evolution of the 
power spectrum. 



3 DATA 

Having developed a suitable formalism for examining the 3- 
D matter power spectrum, we wish to proceed to measure 
this directly from the COMBO-17 survey (Wolf et al 2001). 
This survey currently spans 1 square degree, in four widely 
separated 0.5 x 0.5 square degree fields. One of the fields 
contains the A901/2 supercluster, and is therefore deemed 
highly unrepresentative and inappropriate for inclusion in 
this study (c.f. the significantly anomalous power spectrum 
measured for this field in Brown et al 2003). One further 
field (SGP) does not currently have full redshift informa- 
tion available, so this study focuses upon the two remaining 
fields, centred upon the Chandra Deep Field South (CDFS) 
and an area with no previously known mass concentrations 
(Sll), but which in reality contains a cluster at z = 0.11; 
Brown et al (2003) show that this does not produce a high 
power spectrum estimate in this field. 

All fields in the survey were observed at the MPG/ESO 
2.2m telescope, La Silla, Chile, with the Wide-Field Imager. 
This is a 4 x 2 array of 2048 x 4096 pixel CCDs, with pixel 
scale 0.238 arcseconds (c.f. Gray et al 2002). 

The key characteristic of the COMBO-17 data which 
makes it suitable for our purpose is the existence of an ex- 
tensive catalogue of accurate photometric redshifts for the 
galaxies in each field. Each COMBO-17 field was observed 



in 5 broad-band filters (UBVRI) and 12 medium-band fil- 
ters running from 350nm to 930nm. This permits simul- 
taneous estimates of redshift and SED classification using 
empirically-based spectral templates; this is described in de- 
tail in Wolf et al (2001), who use the same version of the 
catalogue. The resulting accuracy of photometric redshift is 
A z ~ 0.05 for galaxies with z < 1.0. 

The weak shear measurements were made upon the R 
filter images from the survey. This filter was used in the best 
seeing conditions, providing a deep image optimal for weak 
lensing; the combined exposure length was 36120 seconds 
for the CDFS field, and 18100 seconds for the Sll field. 

The reduction process for this R band imaging data is 
described in Gray et al (2002) and Brown et al (2003). For 
the CDFS field, 78 x 8 CCD chip exposures were registered 
by using linear astrometric fits, including a 3a rejection of 
bad pixels and columns; for the Sll field, 44 x 8 exposures 
were registered in this fashion. 

The shear estimation for the galaxies was carried out us- 
ing the imcat package, following the methodology of Kaiser, 
Squires & Broadhurst (1995); see Gray et al (2002) and 
Brown et al (2003) for details of the implementation for 
our dataset. The analysis yields a catalogue of galaxies with 
centroid positions, plus shear estimates which have been cor- 
rected for circularisation of the galaxies by the PSF, and 
anisotropic smearing due to e.g. tracking errors. This cata- 
logue was combined with the photometric redshift catalogue 
for the COMBO-17 survey. 

The resulting catalogue included 52139 galaxies, 23102 
of which had reliable photometric redshifts assigned; the re- 
maining galaxies yielded ambiguous redshift probabilities or 
were fainter than the R — 24 reliability limit for redshifts in 
the survey. This remainder of galaxies was flagged as having 
unreliable redshifts; the faint galaxies were assigned an op- 
tional redshift which we use below when stated. Following 
Brown et al (2003), Section 4, we use a median-magnitude 
median-redshift relation to assign this redshift; the median 
R magnitude of the faint sample is R = 24.9, corresponding 
to a median redshift z ~ 0.95 ± 0.05. 



4 ANALYSIS 

We will now use the models developed above to assess the 
evidence for evolution of the matter distribution from the 
COMBO-17 survey. We begin by making a minimal di- 
rect detection of evolution, by comparing the growth of the 
lensing signal with redshift to that expected with a truly 
unevolving matter distribution. We go on to use our phe- 
nomenological model of an evolving power spectrum to con- 
strain the rate of evolution. Finally, we return to full power 
spectra models including cosmological parameters, in or- 
der to constrain the cosmological constant from the rate of 
growth of the shear signal with redshift. 

4.1 Preliminary detection of evolution 

As a preliminary first step in using 3-D shear data to mea- 
sure evolution, we will demonstrate that the data exclude a 
non-evolving matter distribution, given the present-day nor- 
malisation of the matter power spectrum as,o and a ACDM 
geometry. It is important to realise that this simplest step is 
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already effectively achieved by 2D cosmic shear surveys, if 
the median redshift of the survey is well estimated. That is 
to say, the amplitude of the projected 2D shear field correla- 
tion functions measured with cosmic shear surveys (see e.g. 
van Waerbeke & Mellier 2003, Refregier 2003 for reviews) at 
a median redshift z m ~ 1 is inconsistent with that expected 
for an unevolving power spectrum, both calculated for ex- 
ample with ct 8 ,o = 0.84 ± 0.04 (c.f. Spergel et al 2003). For 
example, at an angular scale of 1', Co ^ (3.2 ± 0.4) x 10~ 4 
for z m — 1 in these surveys, whereas Co is predicted to be 
(5.1±0.8) x 10~ 4 at 1' in a model with unevolving power (i.e. 
with the power expected for ACDM at the present day as in 
section 2.1, but unchanged in the past), cts.o = 0.84 ± 0.04 
and ACDM geometry (Q m = 0.3, Q to t =1,H = 72km s _1 
Mpc -1 ) calculated using equation (19). 

However, we will see that there is considerable power 
in using the redshift information for each galaxy, as we will 
obtain significantly more accurate measurements of the evo- 
lution, and will be able to exclude no-growth models re- 
gardless of present-day power spectrum normalisation. We 
begin by making a direct detection of evolution with the 3-D 
shear field. We can do this by calculating the \ 2 fit of the 
no-evolution model to the data, and comparing with the % 2 
fit for, say, an evolving ACDM model (Q m = 0.3, Qtot = 
l,Ho = 72km s -1 Mpc -1 ), in the first instance both with 
f8,o = 0.84 (c.f. Spergel et al 2003). In the next section 
we will examine the effect of varying this normalisation and 
marginalising over cosmological parameters. 

The sum we require for our \ 2 fit is given by 

X = y^Wij[(7i,»7i,j - Ci(9ij,Zi,Zj)) 2 + 

(72,»72j -C2{6n,Zi,Zj))% (34) 

where 71 and 72 are the two components of shear for a par- 
ticular galaxy, in the frame where the x-axis is the line join- 
ing the two galaxies in question; Wij are a set of weights 
which we will choose presently, i and j are indices for dif- 
ferent galaxies, and Ci and C2 are the expected correlation 
functions between the two galaxies in the model, calculated 
using equations (21) and (22). 

Here we choose a simple weighting Wij = l/e>" 4 , i / 
j, where cr 7 is the variance in one shear component. We 
choose this weighting as the error upon 71^71^ is entirely 
dominated by the random orientation of galaxies; therefore 
we have neglected off-diagonal elements of the covariance 
matrix of errors upon 71,171 j. In addition, we neglect the 
smoothing error associated with our redshift uncertainty of 
Az ~ 0.05; assessment of the impact of this small smoothing 
term is left for future investigation. We are also invoking 
the Central Limit Theorem in order to use \ 2 to estimate 
uncertainties upon parameters, as any one product of galaxy 
shears does not have a Gaussian error. 

In order to remove the impact of intrinsic alignments 
between galaxies, due to physical galaxy alignment instead 
of apparent alignment due to lensing, we reject galaxy pairs 
which are within Az = 0.05 of each other, if their comoving 
separation is less than IMpc. We also remove, with the same 
prescription, near galaxy pairs in the background without 
redshift information, assigning to them the median redshift 
as described in Section 3. 



The weighting equation (34) is 1/cr 4 since, if our shear 
estimator 7; is near Gaussian, we find that the probability 
distribution of 7^ (i.e. two independent Gaussian variables 
multiplied together) is given by 

Prob(j/ = 7i7j -) = -1^K (lL) , (35) 

where Ko is a modified Bessel function of the second kind. 
This probability distribution, Prob(y), has variance <r 4 as 
stated; this is found to be an excellent match to the variance 
of measured 7;7j (both have value 0.31 4 ). 

The no-evolution model we select is as described above, 
i.e. ACDM (On = 0.3,fitot = 1,H = 72km s -1 Mpc -1 ) 
power at the present day with a 8 ,o = 0.84; however, 
this power remains constant at earlier epochs. An evolv- 
ing ACDM geometry is assigned with parameters as above, 
evolving as appropriate to a given epoch; that is to say, the 
lack of evolution is directly in the power, not in the geometry. 
As described in Section 3, we apply the fit to two cases: in 
the first case, we include only galaxies with a reliable photo- 
metric redshift measurement. In the second case, we assign 
a redshift z = 0.95 to objects without measured redshift, 
for both evolving and non-evolving measurements of \ 2 ■ (In 
the next section, we will include the effects of marginalisa- 
tion over the allowed range of median redshift and Hubble 
parameter) . 

If we include the background objects with assigned 
median redshift, we find that the probability ratio be- 
tween the no-evolution model and the ACDM model is 
ex P (Xnoov _ Xacdm)/2 = e -31 ' 1 ; thus the no-evolution 
model is completely excluded given this simple two-model 
comparison. If one only uses objects with known redshifts, 
one finds a probability ratio of 49.4; i.e. the ACDM model 
is approximately 50 times more likely a model than the no- 
evolution model, given our data. 

Thus we see that, using shear estimators together with 
photometric redshifts for each galaxy, we are immediately 
able to set powerful constraints on the evolution of struc- 
ture. But we must now proceed to a fuller treatment for 
measuring the evolution, including the effect of power spec- 
trum normalisation. 



4.2 Measuring the evolution of the power 
spectrum 

In this section, we seek to measure the slope of the mass 
power spectrum in the redshift direction, together with 
the present-day amplitude of the power spectrum. In or- 
der to achieve this, we use the parameterised power spec- 
trum of equation (33), rescaled as Pg(k, w) = Ak a e~ sz , and 
calculate from this the shear correlation functions as de- 
scribed in Section 2.2, using our standard ACDM geome- 
try (Q m = 0.3,fi t ot = I, H = 72 ± 5km s -1 Mpc -1 from 
Freedman et al 2001). We carry out the \ 2 calculation above 
(equation (34)), varying the parameters A (power spectrum 
amplitude at present day) and s (related to the gradient of 
the spectrum as a function of redshift), while fixing a — —1.2 
as discussed in Section 2. We marginalise over Ho = 72±5km 
s -1 and median redshift z ~ 0.95 ± 0.05. 

The resulting constraints on the parameters A and s are 
shown in Figure 1. The panels show the constraints from the 
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Figure 1. Joint constraints on power spectrum amplitude A and 
evolution slope s from COMBO-17 survey (CDFS and Sll fields). 
Top panel: combined constraints from the two fields. Dotted lines 
show the constraints from galaxies with known redshifts only. 
Contours display 1 and 2cr confidence regions. Solid lines show 
constraints including galaxies with unknown rcdshift, assigned 
z=0.95. Middle panel: constraints from the CDFS field only. Bot- 
tom panel: constraints from the Sll field only. 



combined fields, CDFS and Sll respectively; dotted lines 
show the constraints if we only include galaxies with known 
redshifts, while solid lines show constraints including the 
galaxies with unknown redshifts. 

It is easy to understand the basic shape of the confi- 
dence contours; if one has a higher amplitude of power at the 
current epoch, one requires a more dramatic rate of growth 
to account for the shear estimates (c.f. the previous section, 
where we saw that a moderate present-day amplitude abso- 
lutely requires decreasing power as we increase redshift). 

Note that our fit allows very large values of the present- 
day normalisation A and slope s (Figure 1 only shows the 
A — s region of interest; acceptable fit is achieved even for 
A — 100). This is because the lensing predictions from a 
finite field are insensitive to the power and growth of power 
present at low redshift (c.f. discussion concerning Figure 2 
below). Thus for our investigations, we will impose a prior 
A < 10 corresponding to a rate of growth approximately 
three times faster than that expected for ACDM, as dis- 
cussed below. 

We find that s < is excluded at the 1.9a level if we 
include all galaxies, and at the 1.2a level if we only include 
galaxies with known redshifts. This initially appears to be a 
much lower level of confidence to that discussed in the previ- 
ous section; however, here we are still allowing models which 
have very low present-day power spectrum amplitudes. If we 
require A > 3, corresponding to a s > 0.4, we find that s < 
is excluded at the 7.7a level if all galaxies are included, and 
at 2.7a in the known-redshift case; therefore we are again 
finding very strong rejection of the no-evolution model, if we 
require a reasonable present-day power spectrum amplitude. 

Turning now to the actual measurements of A and s 
afforded by our fit, we find that our 2-D contours constrain 
s > 2(la) for CDFS or Sll independently (including galax- 
ies with unknown redshifts) . This is entirely consistent with 
our fiducial ACDM model, which is found to have (by fit- 
ting a tangent to the calculated power spectrum model) an 
initial (z — 0) slope coefficient s = 1.7 for k — IMpc -1 and 
s = 2.0 for k = 50Mpc _1 . 

The amplitudes of the two fields show significant vari- 
ation from field to field, with best fit A = 3.2 and A — 4.5 
at s — 2 for CDFS and Sll respectively (including galaxies 
with unknown redshift); thus when the \ 2 are combined for 
the two fields, the reduced overlap artificially increases the 
acceptable values of s for the combined data. Clearly there- 
fore, we will greatly benefit from applying the methodology 
to many more lines of sight. 

A possible extension of this approach is to directly 
couch A in terms of as; however, since our phenomenological 
power spectrum model has a different k and z dependence 
to a full cosmologically parameterised power spectrum, val- 
ues of as calculated for our model do not have the conven- 
tional calibration. Instead, we will carry out the equivalent 
measurement of the normalisation of the power spectrum in 
terms of as in section 4.3. 

One should note the extraordinary improvement in ac- 
curacy arising from using our whole sample, as opposed to 
only those galaxies with known redshifts. The combined con- 
straint upon A and s has improved by a factor of > 2 for 
each field. This acts as a driving concern for future surveys; 
we would clearly wish to have a sample with known red- 
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shifts to R ~ 25, which poses a serious challenge to current 
photometric redshift techniques. 

We can use the information we have acquired on A and 
s to plot constraints upon the power spectrum growth itself. 
Since we can describe Pg = Ae~ 3Z k a , and since we have 
measured the probability distribution of A and s, Prob( J 4, s), 
we can calculate the probability of any particular value of 
Pg at a given z, Prob(P^). 

We calculate this by examining the cumulative proba- 
bility distribution F for Pg, i.e. the probability that Pg is 
smaller than a certain value Pgo- We see that 

F= Prob(Pg)dPg= / dA dsProb(A,s).(36) 

Jo Jo Jllog(^) 

Now Prob(P a ) = dF/dPg, so 

ProMP^^-^Prob^-ilog^)). (37) 

We can calculate this directly from our (A, s) proba- 
bility contours. We plot the resulting constraints upon the 
power spectrum growth in Figure 2. Here, we have used the 
probabilities including the galaxies without known redshifts, 
as we found above that the constraints from the redshift 
sample are weak. Also plotted is the growth prediction for 
ACDM (Qa = 0.7, On = 0.3) with as = 0.7 (c.f. section 4.3); 
we see that the A 2 constraints are consistent (within 2a) 
with this model, for < z < 1. Agreement is better between 
the model and an individual field, rather than the fields com- 
bined, as discussed above. Note that a k — 14Mpc _1 has 
been chosen for this plot; the phenomenological model used 
presents us with an average, fe-independent redshift evolu- 
tion, as the model has k and z as independent variables. 

The growth of the power spectrum shown on Figure 2 is 
skewed towards achieving high amplitude at the present day, 
with corresponding large s, simply because our small dataset 
cannot currently differentiate well between slow-growth low- 
amplitude and rapid-growth high-amplitude models. Since 
there is a large area in parameter space representing the lat- 
ter, the overall constraints favour these while not excluding 
slower, ACDM growth. Larger surveys will be required in or- 
der to improve the measurement of matter power spectrum 
growth, by measuring the growth of the shear correlation 
function more accurately as a function of redshift and thus 
restricting the A — s degeneracy. With these larger surveys, 
the procedure we have outlined will allow detailed study of 
the collapse of structures in the Universe. 

In this section, we have measured the growth rate 
and amplitude of the matter power spectrum using a phe- 
nomenological model. However, this is only one approach to 
the extraction of useful information from a 3-D shear field; 
instead of fitting a phenomenological model as above, we can 
directly fit the shear predictions for a range of cosmological 
parameters. We will carry out this analysis in the following 
section. 

4.3 Constraints on the cosmological constant 

Instead of a direct parameterisation of the slope of the power 
spectrum with redshift, we can characterise the 3-D shear 
data by the constraints which they provide on cosmolog- 
ical parameters. With the current small dataset, it is not 
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Figure 2. Constraints on power spectrum as a function of red- 
shift z, displayed for k = 14Mpc _1 . The contours show 1 and 2a 
confidence for the behaviour of the power spectrum at a given 
redshift. Also plotted (dashed line) is the prediction for ACDM 
with as = 0.7. Top panel: results for CDFS and Sll combined; 
middle panel: results for CDFS only; bottom panel: results for 
Sll only. 
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possible to obtain cosmologically interesting constraints on 
many parameters at once. Instead, here we will use priors 
from WMAP for most parameters, and examine the result- 
ing constraints on the normalisation erg and cosmological 
constant A from the COMBO-17 data. Future work with 
large 3-D shear datasets will be able to obtain estimations 
of full parameter sets from shear alone. 

4.3.1 N on- flat Models 

We choose in the first instance to use the Hubble Key 
Project and WMAP measurements of parameters Ho = 
72 ± 5km s _1 Mpc -1 and Q b = 0.047 ± 0.006 (Freedman 
et al 2001, Spergel et al 2003), and examine the slice in 
parameter space at Q m = 0.3 allowing Qa and erg to vary 
(i.e. we will allow non-flat models, Qa + Qm 7^ 1). This will 
enable us to examine the constraint from our 3-D lensing 
method on Qa if Q m is known. We calculate the 3-D shear 
correlation functions for this set of parameters as described 
in Section 2. We then make \ 2 fits to the data as above, 
varying erg in steps of 0.05 and Oa in steps of 0.1. 

The resulting weak constraints on these parameters, 
while either excluding or including the unknown redshift 
sample fixed at 2 = 0.95 (with ±0.05 marginalisation to ac- 
count for the uncertainty in median redshift, plus marginali- 
sation over Ho), are shown in Figure 3. In the case where we 
use only galaxies with known redshift, we obtain only very 
broad constraints on erg and Qa- We find an upper limit on 
Qa of Qa < 4.0 — 4.2(78 at the la level, or Qa < 5.15 — 4.8crg 
at the 2a level. These results also constrain as to have a 
low normalisation, as < 0.70(0.93) at the l(2)er confidence 
level for ft a = 0.7. This is consistent with the as normal- 
isation found for the COMBO-17 dataset for a 2D cosmic 
shear analysis (Brown et al 2003). Perhaps this is not sur- 
prising for the two fields of interest here, as both are rather 
empty fields containing little evidence of substantial large- 
scale structure, other than the cluster in one chip of the Sll 
field. Therefore further fields will require analysis before one 
can obtain a definitive value for as from this method. 

If we include galaxies with unknown redshifts, assigned 
a redshift 0.95 ±0.05, our constraints remain weak but show 
that the Qa estimate does vary weakly with as; the best 
fit joint constraint line is Qa = 5.6 — 7.7(crg ± 0.07). This 
behaviour is as expected; the weakness of the dependence 
shows that ACDM and OCDM models have quite similar 
evolution over < 2 < 1. As seen above, these results give 
a low estimate of as on our fields; note that we already 
expected these fields to give a low estimate for the spatial 
fluctuation of dark matter. 



4.3.2 Flat Models 

An alternative means of deriving useful constraints from this 
small dataset is to fix Qa + fl m =l, i.e. we can impose the 
flatness condition measured by WMAP (Qtot = 1-02 ± .02 
with Hubble key project Ho prior, Freedman et al 2001, 
Spergel et al 2003). With this condition, we can again make 
shear correlation function predictions, then obtain \ 2 fits to 
the data for various values of as and Qa, incrementing in 
units of 0.05 and 0.05 respectively. We marginalise over Ho 
and ^median as before. 




Figure 3. Constraints on Qa and as from the COMBO-17 sur- 
vey, with fi m fixed to be 0.3 and curvature allowed to vary ar- 
bitrarily. Dotted lines (and lightest two greyscales): constraints 
solely from sample with known redshifts. Solid lines (and darker 
two greyscales) : constraints including galaxies with unknown red- 
shift. 



The results are shown in Figure 4, for the cases where 
we exclude and include the unknown redshift sample fixed 
at 2 = 0.95 ± 0.05. Here we find that the constraints on Qa 
or Qm are substantially improved over the case where we do 
not prescribe flatness. 

In the case where we exclude the unknown redshift 
galaxies, we find a weak constraint Qa > 1 — 0.16a g _1 ' 95 (la) 
or Qa > 1 — 0.25ag 17e (2a). If we then force a prior 
as = 0.84 ± 0.04, we obtain Q A = 0.91±q°| 7 at the 2a level, 
consistent with standard ACDM as defined above. 

Again, these results offer a low normalisation of the 
matter power spectrum; for Qa = 0.7, we require as < 0.72 
at the la level. As discussed in the previous section, this re- 
sult is dominated by the fact that our two fields are devoid 
of significant large-scale structure; further fields will require 
measurement before cosmological implications can be con- 
clusively drawn. 

If we now include galaxies with unknown redshift, we 
find a best-fit constraint Qa = 1 — 0.15(crg ± 0.04)~ 1,5 , erg < 
0.72(ler). In the case where we impose a prior as = 0.84 ± 
0.04, we obtain Q A = 0.83lg;??(2cr). This, and the constraint 
directly from our contours Qa < 0.76(la) is in accord with 
the concordance ACDM model; only the power spectrum 
normalisation is outside the range expected. 

Despite the limitations due to the size of our dataset, 
these constraints are already impressive, and demonstrate 
the power of using 3-D information. Figure 4 compares the 
accuracies of the 2D analysis of Brown et al (2003, dashed 
lines) for only the CDFS and Sll fields with the current 3-D 
analysis (solid lines) ; clearly one gains very significantly from 
the inclusion of the redshift information, by a factor of > 2 
everywhere in uncertainty (beyond as > 0.4). This holds out 
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Figure 4. Constraints on Qa and eg, with Q\ + f2 m fixed to 
be 1. Dotted line (and lighter two greyscales): constraints solely 
from sample with known rcdshifts. Solid line (and darker two 
greyscales): constraints including faint sample. Dashed lines show 
the constraints from our 2D lensing analysis in Brown et al (2003). 

the promise of very precise measurements of cosmological 
parameters with future 3-D lensing surveys. 



5 CONCLUSIONS 

In this paper, we have explored the evolution of large-scale 
structure for redshifts z < 1. This has been achieved using 
weak gravitational lensing together with photometric red- 
shifts for a set of galaxies, both measured upon two fields 
from the COMBO-17 survey (Wolf et al 2001). 

We have described the construction of theoretical mod- 
els of evolving matter power spectra, both in terms of phe- 
nomenology of the growth of structure, and in terms of char- 
acterisation of this growth with cosmological parameters. 
We have then shown how to calculate the shear power spec- 
trum from this matter power spectrum, and vice versa. In 
particular, we have discussed the importance of the cross 
power spectrum for shears between redshift shells, and have 
related phenomenological models for evolution between the 
matter and shear cross power spectra. 

In order to practically implement this theoretical frame- 
work, we have described the usefulness of the COMBO-17 
survey for 3-D shear analysis, due to its accurate photomet- 
ric redshifts for galaxies in the survey (with Az ~ 0.05 for 
< z < 1.0) and its well-studied shear catalogue. 

We proceeded to apply a least squares estimator anal- 
ysis to the data, finding that best-fit matter power spec- 
tra included significant evolution. In particular, we found 
that the growth rate expected for ACDM was approximately 
50 times more likely than a zero evolution model, in the 
case where we ignore galaxies with unknown redshifts. How- 
ever, we found that the best-fit rate of growth understand- 
ably depended sensitively upon the present-day amplitude of 



the power spectrum. We calculated constraints upon A 2 (z), 
finding a reasonable match with that predicted by ACDM 
for < z < 1, but allowing larger present-day amplitude of 
the power spectrum with more dramatic growth rates. 

As an alternative route to extracting information from 
the 3-D shear field, we obtained least squares estimates 
of the cosmological constant and the normalisation of the 
power spectrum. We find low values of as in this small area 
of sky, with best-fit constraints Q,a = 5.6 — 7.7(o"8±0.07) and 
Oa = 1 - 0.15(ct8 ± 0.04)- 15 ,fT 8 < 0.72 if we use WMAP 
priors on Q m or ilk respectively. We find that the latter 
constraint is a factor of 2 improvement upon the constraint 
obtained by 2D lensing for the same area of sky. 

These results demonstrate that 3-D statistical lensing is 
a useful means of examining dark matter evolution. If one is 
interested in measuring only cosmological parameters, this 
method is a complementary alternative to measuring 2-D 
higher-order statistics of the shear field; the question of the 
comparative (and combined) merits of the two approaches is 
the subject of ongoing research (c.f. Benabed & van Waer- 
beke 2003, Refregier et al 2003). However, if one seeks a 
direct measurement of the dark matter evolution itself, 3-D 
methods such as the one explored here are essential. With 
large multi-colour imaging surveys in the next few years such 
as that of SNAP (Rhodes et al 2004), the potential of this 
approach will be fully realised, affording measurements of 
the cosmological constant and dark energy equation of state 
parameters to unprecedented accuracy, together with pre- 
cise reconstructions of the growth of the power spectrum as 
a function of redshift. 
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